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transform Program 
for 
colour 
coordinate 

finclude <stdio.h> 

void 

oainO 

{ 

float 

dWCaj.alSlCSl.bCSl.cCal.ddCal^r^Kl.kkCaltal.ss.sss; 
int i,j»k; 

/* input z & 7 of RGBW */ 

printfCIHPDT RGB and White¥n'); 

printfCEx Ry Gx Gy Bx By Wx tfy¥n*); 

scanfCXf Xf # %f Xf Xf Xf XT, td[0][0],*d[0][l] 

, 4d[l][0],»d[l][l] 

, &d[2][0],&d[2][l] 

, ftd[3][0],«[3][l]); 

/* 



d[0][0] 




0.67; 




d[0][l] 




0.33; 




d(l][0] 




0.21, 




dtl][l] 




0.71, 




d[2][0] 




t).14, 




d[2][l] 




0.08, 




d[3][0] 




0.31, 




d[3][l] 




0.316; 



*/ 
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/* calculate z fron z & y */ 
for(i = 0; i < 4; i++){ 

if((d[i][0] + d[i][l]) > 1.0){ 
dli][2] = 0.0; 

} 

d[i][2] = 1.0 - d[i][0] - d[i][l]; 

} 

printfCHATBIHn"); 
for(i = 0; i < 3; i++M 

printf('¥f ); 

for( j = 0; j < 3; j++){ 

printt("X5.3f¥t" ) d[i][j]); 

} 

printf("¥n»); 

} 
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/* caluculate matrix */ 
{ 

int il, iZ, jl, j2; 
for(i = 0; i < 3; i++){ 

11 = i + 1; 

12 = i + 2; 

if (il > 2) il = 0; 
if (i2 > 2) i2 = i2 - 3; 
for(j = 0; j < 3; 
jl = J + l; 
j2 = j + 2; 
if (jl > 2) jl = 0'; 
if (j2 > 2) j2 = j2 - 3; 

a[i][j] = d[il][jl]*d[i2][j2] - d[il][j2]*d[i2][jl]; 

} 

} 

} 

/* calculate of BlfflBO */ 
for(i = 0;i < 3; i++){ 

b[i] = 0; ■ 

for(j = 0; j <3; j++){ 
b[i] = a[i][j] » d[3][j] + b[i); 

} 

} 
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/* MATRIX */ 
for(i = 0; i < 3; i++){ 
for(j = 0; j < 3; j++){ 
a[i][j] = a[i][J] / h[i]; 
r[i][j] = a[i][j]; 
if(i = j){ 

kk[i][j] = 1.0; 
} else { 

kk[i][j] = 0.0; 

} 

} 

} 

/* IHVEfiSE MATRIX */ 
for(i =0; i < 3; i++){ 
for(j = 0; j < 3; 

ddCj] = a[j][i]; 

a[j][i] = 0.0; 

} 

a[i][i] = 1.0; 
for(j = 0; j < 3; j++){ 
c[j] = a[i][j] / dd[i]; 

} 

for(j = 0; j < 3; J++M 
for(k = 0; k < 3; k++){ 
a[j][k] = a[j][k] - c[k]*dd[j]; 

■ } 
} 

for(j = 0; j < 3; j++H 
a[i][j] = c[j]; 

} 
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/* SEIKIKA */ 

ss = a[l][0] + a[l][l] + a[l][2]; 
sss = r[l][0] + r[l][l] + r[l][2]; 
for(i = 0; i < 3; i++){ 

for(j = 0; j < 3; j-H>{ 

a[i][j] = a[i][j] / ss; 

r[i][j] = r[i][j] / sss; 

} 

} 
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/* result */ 

printf (" original data¥iT ) ; 
for(i = 0; i < 4; i++){ 

printf("¥t B ); 

for( j = 0; j < 3; j++M 

printf ("J7.5f ",d[i][j]); 

} 

printfCMta"); 

} 

printfrHAraiHB"); 
for(i = 0; i < 3; i++){ 

printf("¥t"); 

for( j = 0; j < 3; j++){ 

printf (-%7.5f -,r[i][j]); 

} 

printf W); 

> 

printf (" INVEECE HATRIZSn" ) ; 
ford = 0; i < 3; i++){ 

printf(-¥t-); 

for( j = 0; j < 3; j++){ 

printf("X7.Sf -,a[i][j))-; 

} 

printf("¥n'*); 

} 

for(i = 0; i < 3; i++){ 
for(j = 0; j < 3; 
tt[i][j] = a[i][0]*r[0j[j] + a[i][l]*r[l][j] + a[i][2]*r[2][j]; 

} 

} 

printf("KAKEZAN¥n**); 
for(i = 0; i < 3; i++){ 

printf C¥t"); 

for( j = 0; j < 3; 

printf("%7.5f -,kk[i][j]) : 

} 

printfCTfn"); 

} 
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